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Factorization in B — > Kir£ + £ decays 
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We derive factorization relations for the transverse helicity amplitudes in B — > Kir£ + £~ at leading 
order in A/mj, in the kinematical region with an energetic kaon and a soft pion. We identify and 
compute a new contribution of leading order in A/rrib to the B — > K-r£ + £~ amplitude which is not 
present in the one-body decay B — > K*£ + £~ . As an application we study the forward-backward 
asymmetry (FBA) of the lepton momentum angular distribution in B — > Ktt£ + £~ decays away 
from the K* resonance. The FBA in these decays has a zero at q$ = q^MKir), which can be 
used, in principle, for determining the Wilson coefficients CV,9 and testing the Standard Model. We 
point out that the slope of the qo(Mx w ) curve contains the same information about the Wilson 
coefficients as the location of the zero, but is less sensitive to unknown nonperturbative dynamics. 
We estimate the location of the zero at leading order in factorization, and using a resonant model 
for the B — > Kir£ + £~ nonfactorizable amplitude. 

PACS numbers: 12.39.Fe, 14.20.-c, 13.60.-r 



I. INTRODUCTION 



The rare electroweak penguin decays fo — >■ 57 and 
b — > s£ + £r arc sensitive probes of the flavor structure 
of the Standard Model, and provide a promising testin 
ground for the study of new physics effects (sec Rcf. 
for a recent review of the experimental situation). Sev- 
eral clean tests have been proposed in these decays, which 
are sensitive to the chiral structure of the quark couplings 
in the Standard Model. Examples of such tests involve 
measuring the photon polarization in b — > 37 and the 
zero of the forward-backward asymmetry in b — > s£ + £~ 







Our understanding of these decays has advanced con- 
siderably over the past few years, through the derivation 
of factorization relations for exclusive B — > K^*H + £~ 
and B — > K*j decays at large recoil. First derived 
at lowest order in perturbation theory 0, , these fac- 
torization theorems were proved to all orders in a s 
j3H|I|i^EiEll2i| using the soft-collinear effective 

In this paper we introduce a new factorization relation 
for the multibody rare decays B — ► Kir£ + £~ in the kine- 
matical region with a soft pion and an energetic kaon, 
at leading order in A/mi,. The schematic form of the 
factorization relation is given below in Eq. (|40|l . This 
extends the application of factorization to final states in 
B — > X s £ + l~ containing a few hadrons, with small total 
invariant mass. 

A particularly clean test for new physics effects in these 
decays is based on the forward-backward asymmetry of 
the (charged) lepton momentum in B — > K*£ + £~ with 
respect to the decay axis q = pi+ + p%- . This is defined 



as 



A FB {q 2 ) 



dT{q 2 )/dq 2 



dcosf 



dq 2 d cos 9+ 

dT(q 2 ,9 + ) ' 
dq 2 d cos 9 + 



(1) 



where 9+ is the angle between p e + and q in the rest frame 
of the lepton pair. 

As pointed out in HB, ^ uc t° certain form factor 
relations at large recoil 0, H| , this asymmetry has a zero 
at (?g which depends mostly on Wilson coefficients in the 
weak Hamiltonian with little hadronic uncertainty. The 
position of the zero was computed in 0,0,0, using the 
complete leading order factorization formula. The most 
updated result, including isospin violation effects, is [Tfij 



% 



Re(Cf(q 2 )) 



Rc {Cf{q 2 )) 

4.15 ±0.27 GcV 2 (B+ 
+0 il GeV 2 (B° 



(l + r(a s ))+c5 fact (2) 



4.36 



K*+) 
K*°) 



Here r ~ a s (rrib) is a radiative correction and Jf ac t de- 
notes factorizable corrections which break the form factor 
relations. Precise measurements of the position of the 
zero <7q can give direct information about new physics 
effects through the values of the Wilson coefficient C| ff 
(with Cj S determined from B ~ > X s -f decays). 

The branching ratios of the B — > K*£ + £~ exclusive 
modes have been measured with the results 



B(B -> K*£+r) = 

(7.8^7 ±1.2) x 1(T 7 
(16.5 ±2.3 ±0.9 ±0.4) 



(3) 



x 10" 



(BABAR) 
(BELLE) 



Differential distributions of the g 2 -spectrum are also 
available, as binned branching ratios. First measure- 
ments of the forward-backward asymmetry Apsiq 2 ) have 
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been presented by the BELLE Collaboration [5j , but due 
to large errors the position and even existence of a zero 
are still inconclusive. 

In practice, the K* is always observed through its 
strong decay products K* — > Kir. We point out that 
this has several interesting implications. The multibody 
factorization relation proven here contains a new factoriz- 
able contribution to the decay amplitude of leading order 
in A/mi, which is not present in the B — > K*£ + £~ factor- 
ization relation. This introduces a shift in the position 
of the zero of the FBA in this region. As an applica- 
tion of the new factorization relations, we compute the 
correction to the position of the zero arising from this 
effect. 

Another novel effect is the existence of a zero of the 
FBA also for a nonresonant Ktt pair, which occurs at a 
certain dilepton invariant mass q^Mjen) depending on 
the hadronic invariant mass Mktt- In principle, this ex- 
tends the applicability of the SM test using the zero of 
the FB asymmetry also to nonresonant B — ► Kir£ + £~ 
decays. In practice however, the calculation of the posi- 
tion of the zero is complicated by the appearance of ad- 
ditional nonpcrturbative contributions to the amplitude. 
We estimate the unknown nonfactorizable amplitude in 
B — > Kir£ + £~ in terms of a K* resonant model. 

We propose an alternative test of the SM using the 
slope of the q^(M^) curve, which can be shown to con- 
tain the same information about the Wilson coefficients 
as the location of the zero itself. In contrast to the ab- 
solute position of the zero, which depends on less well 
known hadronic parameters, the slope of the zero curve 
can be shown to be less sensitive to such effects. 

In Sec. II we introduce the SCET formalism and write 
down the effective Lagrangian for the rare B — > X s e + e~ 
decay. Sec. Ill presents the factorization relations for 
the B — > Kir£ + £~ helicity amplitudes in the kinematical 
region with a soft pion and a hard kaon. Sec. IV lists the 
expressions for distributions in these decays, and gives 
a qualitative discussion of the zero of the FBA in the 
nonresonant region. Sec. V contains a numerical analysis 
of the asymmetry, and finally Sec. VI summarizes our 
results. 



II. SCET FORMALISM 

In the Standard Model the AS = 1 rare B -> X s £ + t 
decays are mediated by the weak Hamiltonian 



H 



G 
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We use everywhere in this paper the operator basis for 
Oi_6 defined in Rcf. [24|]. Smaller contributions to the 
amplitude arise from T-products of the operators in 
Eq. J2J with the electromagnetic current. 

We choose the kinematics of the decay such that the 
total dilepton momentum <j M = + Pi-)^ points 

along the — e$ direction, and has the components q = 
(q , 0, 0, — \q |), expressed in usual four-dimensional coor- 
dinates a/ 1 = (a , a). The hadronic system moves in the 
opposite direction +S3 in the B rest frame. We define the 
light-cone unit vectors n/ 1 = (1,0, 0, 1), ft/ 1 = (1, 0, 0, -1). 
They can be used to project any vector a M onto light- 
cone directions, according to a+ = n^a^ and a_ = n^a^ . 
Finally, we introduce a basis of orthogonal unit vectors 
e± = ^(0, 1, Ti, 0), £ = 0, 0, q ). 

We will be interested in the kinematical region with 



q < mi 



for which the hadronic system has a large 



light-cone momentum component along n. This defines 
the hard scale Q = n- px ~ rrib ^> A, with A ~ 500 MeV 
the typical scale of the strong interactions. 

The effective Hamiltonian Eq. is matched in the 
SCET/ onto 



G 



V2 



K s} - {{h^)J v + {HiM} (6) 



where the currents J v A have each the general form 



...(<) 



(7) 



+ [6g(^K + &S( Wj -K] J (10) K) 

with i = V,A, and integration over u>j = (u>i,L02) is im- 
plicit on the right-hand side. This expansion contains the 
most general operators up to order O(A), with A 2 = A/Q. 
The subleading operators are defined as 



t(il,lr)/, , , , \ _ - r( 1L - 1 
• J n '(^1,^2) - Sn,D,lL 



R) 



1 



■AgBlr. 



b r . 



Wi,w 2 J = q n ,ui 



1 



PL br, 



with 



r r (lL) p(lfl)-. 



n ■ V 



(8) 



The 

collincar gauge invariant fields are defined as q n = W'£ n , 
igBfj, = W^[n ■ iD c , iD^W, with £ n the collinear quark 
field and IF = cxp[—g(n-A ntq )/(n-q)] a Wilson line of the 
collinear gluon field. We use throughout the notations of 
Ref. |l9j | with n ■ v = n ■ v = 1. 

The Wilson coefficients of the leading order SCET op- 

(4) 

erators appearing in the matching of Jy A are |l6l | 



with Ao 



VtbVf*. The dominant contributions come 



from the radiative penguin O7 = j% so 'pvPnF^b and 
the two operators containing the lepton fields £ = e, fi 

9 = -(s^P L b)(£r£) , do = -(s^P L b)(£^j 5 £) .(5) 

7T 7T 



x(1 _^F^ (W)/i))+0(a 2 (Q)) 



(9) 



c[ A \^n) = C w (l + 0(a 2 s (Q)) (10) 
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The effective Wilson coefficients Cj g include the contri- 
butions of the operators Oi_6 and O s . For convenience 
they are listed in the Appendix, together with the func- 
tions f v (u),fj,) containing the 0(a s (Q)) contribution to 
the Wilson coefficient of the vector current in SCET, and 
t(uj, fi) giving the additional contribution from the tensor 
current. 

The Wilson coefficients of the O(A) SCETi operators 
are given at leading order in a s (Q) by 



2n • q 



(wi,cj 2 ) = [2m b (n) 



dv) 

7 IR 



C L - 



8mb 
-Cf + C 



e u C 2 t±(x,m c ) 
1 



2n ■ q 



(wi,w 2 ) 







0^o{wi,U2) = C10 — 



(11) 

(12) 
(13) 



with uji = xuj ,u>2 = —xlo, and x = 1 — x. We neglect 
here smaller contributions from the operators O3-6 and 
the gluon penguin 0%, which will be retained only in 
the leading order SCET Wilson coefficients c^' \u>, u). 
The complete expression can be extracted from Ref. [3J. 
The Wilson coefficients Ci are defined in the Appendix. 
The function t±(x, m c ) appears in matching from graphs 
with both the photon and the transverse collinear gluon 
emitted from the charm loop and is given in Eq. I|A.8|) 
of the Appendix. 

The coupling of the virtual photon 7* — > to the 

light quarks can also occur through diagrams with inter- 
mediate hard-collinear quarks propagating along the pho- 
ton momentum |26j| . (Such a description is appropriate 
only for a range of the dilepton invariant mass q 2 < 1.5 
GeV 2 which can be considered hard-collinear, see below.) 
Such contributions are mediated by new terms in the 
SCETi effective Lagrangian, which in the SM contains 
only one operator at leading order 



4Gf a ( s) 



n/2 



(14) 



n,uj 3 ifiP L b v )(q 

1 1 . UJ 1 



The Wilson coefficient is given by 



>(«) 



A 



0) 



(c 2 + ^)v-(c 4 + |; 



(15) 



-0(a s (Q)) 



with the soft spectator quark in the B meson [TtI Il8[ . 
Working to leading order in SCET n, t his matching con- 
tains two types of operators [1 ITllia EH. 



J? 



4 4) (c)^]o; f 



(16) 



The first type are the so-called 'nonfactorizable' opera- 
tors, denoted here as 0£ { ,0' n{ . They are defined such 
that they include the contributions of the leading op- 
erators in the SCETi Lagrangian. The second type of 
operators are the so-called factorizable and spectator in- 
teraction operators, denoted as and J^ p , respectively. 
Although their form is similar, they arise in matching 
from different operators in SCETi, as follows. The op- 
erators Jf are obtained from the O(A) operators in the 
SCETi current. The spectator operators contribute only 
to Jy in Eq. ©, and arise from the SCETi weak nonlep- 
tonic effective Hamiltonian Eq. i|14|) ■ The ellipses denote 
terms suppressed by powers of A/m&. 

The matrix elements of the nonfactorizable operators 
in Eq. corresponding to a B — > M„ transition, arc 
parameterized in terms of soft form factors. We define 
them as @ 



(M n ( PM 
(M n (p M 



\0' nl \B) = 



2EmCo(E m ,h) (17) 
)=2E u a(E M ) (18) 



where in the first matrix element M n is a pseudoscalar 
meson, and in the second M n is a transversely polarized 
vector meson. 

The factorizable operators in Eq. (|16f) arc nonlocal 
soft-collinear four-quark operators. As mentioned above, 
they are of two types, denoted as factorizable- type (f), 
and spectator- type (sp). The J? { operators have the 
generic form 

Jft~ J dxdzdk + b {i) (z)J(x,z,k+) (19) 
x(q k+ T s b v )(q ntUJ1 Tcqn,u 2 ) 

where b(z) are SCETi Wilson coefficients, and J are jet 
functions. They are given in explicit form in Eq. (|29fl be- 
low. We use a momentum space notation for the nonlocal 
soft operator, defined by 



/ £e-***+9(A£)y(A,0)6j(0) (20) 



where we neglect again smaller contributions propor- 
tional to C§. 

For application to exclusive B — + M form factors, with 
M = 7T , p , ■ ■ ■ a light meson, the SCETi effective La- 
grangian in Eq. © has to be matched onto SCETn op- 
erators. This requires taking into account the interaction 



with F(A,0) = Pexp(ig J dan ^ (an/ '2)) a soft Wilson 
line along the direction n^. 

There are two jet functions, defined as Wilson coef- 
ficients appearing in the matching of T-products of the 
SCETi currents Eq. JSJ with the ultrasoft-collinear sub- 
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leading Lagrangian C$ , onto SCET n EHHIp 
Tiqn,»j9B^ 2 ] ia (0)[ig?Un]i b = (y) = 
zS ab 5(y + )S^(v^- C dx [ ^± e ^ + v-/2 

U Jn J 47T 



(21) 



x{J\\(x,z,k + )[(yt'y<ty i lq n ^ LlJ ^q n _ m J\ 

+J±(x, z, k^){i^^) 3% [q n>xu -^q n - muj ]} 

with u>i = zlu and lu = lo\ — dj 2 - The jet functions are 
generated by physics at the hard-collinear scale ^ 2 C ~ QA, 
and have perturbative expansions in a s ( pL r ). At lowest 
order in a s (ii c ) they are given by [l8Lll9j | 

J ul (x,z,k + ) = —^5{x-z)— 22 

11 N c xk + 

Finally, the spectator-type operators have the form 



jf 1 



dzb S p(z) / dk-J sp {k- ) 

o J n-q 

-PLq n -zuj) 



(23) 



The jet function J sp (fc_) is the same as the jet function 
appearing in the factorization relation for B — > ^tv. It 
can be extracted from the results of Ref. [2^| and is given 
at one-loop order by 



Jsp{k—) 



1 



k- + ie 



1 



4tt 



■{V 



1 



with L = log[(— n • gfc- — ie)/fj, 2 ]. 

The matrix elements of the factorizable and specta- 
tor operators in the B — > M n transition at large recoil 
are computed as convolutions of the product of collinear 
and soft matrix elements. Adding also the nonfactoriz- 
able contribution, the generic form of the factorization 



relation for the hadronic matrix element for B 
leptons is written as 01 (with i = V, A) 



{M n \J x \B) = cf ] 2EuC 



M„ + 
(25) 



dxdzdk + b^ (z)Jj(x, z, k + ) 
/ dxb sp (x)(j) M (x) [ dk-J sp (k-)(0\q k _r s b v \B). 



The nonperturbativc soft and collinear matrix ele- 
ments appearing in this relation are given by the B-meson 
and light meson light-cone wave functions, respectively. 
We list here their expressions, adopting the following 
phase conventions for the meson states 

(tt + , 7T°, 7r~) = (ud, —=(uu — dd), du) (26) 



V2 



(B~,B°) = (bu,bd). 



(a) 
b 



n 



q n 



i 
i 

A t 

! ! k+ q b 

=+333333333^ *— = 



(b) \s n q„/ 

\ / 
\ / 
\ / 

=• — >- 



q(k) 



FIG. 1: Leading order SCETi graphs contributing to the 
factorizable amplitude for B — > K n n + leptons. a) the T- 
product {J VA ,iC^}, where the filled circle represents the 
SCET current J v A ; b) spectator-type contribution, with the 
virtual photon attaching to the light current quark; the inter- 
mediate quark propagator can be either hard-collinear along 
the photon direction n M , or hard, depending on the ratio 
n-qA/q 2 being larger or less than 1, respectively. After match- 
ing onto SCEFn these graphs contribute to the operators Ji,f 
and J sp , respectively. 



The light mesons' light cone wave functions are given by 



(KnipK^Sn^yLPLqn^lO) = -fKn-p K <t>K{x) (27) 



and the B meson light-cone wave function is defined as 



2] 



}\ul + bi\B-(v)) = 



(28) 



j i 



The matrix element of the q k _ b v appearing in the last 
term of Eq. I|25|l can be obtained from this by the sub- 



stitution <-> n M . 



A. Factorization in multibody B decays 



We consider here the application of the SCET formal- 
ism to rare B — > K n Tr£ + £~ decays into final states con- 
taining one energetic hadron K n and a soft hadron 7r. 
The heavy-light currents in Eq. © contribute to such 
processes again through T-products with the ultrasoft- 
collincar sublcading Lagrangian |l7j . Typical diagrams 
in SCETi contributing to these T-products are shown in 
Fig.IU 

After integrating out the modes with virtuality p\ c ~ 
AQ connected with the hard-collinear degrees of free- 
dom, these T-products are matched onto SCETn fl8j . 
Performing a Fierz transformation of the four-quark op- 
erators, one finds the following result for the factorizable 
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operators 

J Uact = -^ / dxdzdk + b^ L (z)J ± (x,z,k + ) (29) 

— — / dxdzdfc + 6j 4 ^(z) Jii (a;, z, fc + ) 
2w 7 

-— [ dxdzdfc + [fo^(z)w M + 6^(z)?i^)] 

X J||(x, z, k + ){q k+ 4P L b v ){s ntUJl -P L q ntUJ2 ) 

The coefBcicnts 6j(z) are related to the Wilson coeffi- 
cients &j(wi, ui 2 ) of Eqs. (fTTf) as bj(z) = bj((l — z)u, zuj); 
the labels of the collinear fields are parameterized as 
loi = xu/, u>2 = — w(l — x), lo = u>i — u>2 — n ■ pm, with 
Pm the momentum of the collinear meson M n produced 
by the collinear part of the operator. 

Finally, the spectator-type factorizable operators arise 
from diagrams where the photon attaches to the light 
current quark Fig. \I]p. Although the spectator quark is 
not involved in these contributions, we will continue to 
use the same terminology as in the B — > M n case, due to 
the similarity of the corresponding operators. 

The effective theory treatment of these contributions 
depends on the relative size of the virtuality of the photon 
q 2 and the typical hard-collinear scale rribA. These two 
scales correspond to the two terms k_ + q 2 /n ■ q + ie in 
the propagator of the intermediate quark in Fig. \T]p. 

Several approaches are used in the literature to deal 
with these contributions, which we briefly review in the 
following. One possible approach, used in in QCD 

factorization, is to keep both terms in the propagator, 
and not expand in their ratio. From the point of view of 
the effective theory, this approach is equivalent to treat- 
ing the photon as a hard-collinear mode moving along 
the photon direction (2(|. This approach is certainly 
appropriate for real photons, and for hard-collinear pho- 
tons q 2 ~ 1.5 GeV 2 . It is not clear whether it can be also 
applied to photons with q 2 ~ 4 GeV 2 , as is the case here. 

In this approach the spectator-type effective Hamil- 
tonian Eq. ijTH contributes to exclusive decays through 
T-ordered products with the leading order SCETi La- 
grangian describing photon-quark couplings |26l l29l | . Af- 
ter matching onto SCETn, these T-products are matched 
onto one single operator, which can be written as an ad- 
dition to Jy, and is given by 

1 q ^ d ,s Jo P J n-q 
i 

) (30) 

For consistency with the other factorizable operators in- 
cluded, we work to tree level in a s {Q) 1 but keep terms 



of 0(a s ([j, c )) in the matrix elements of the factorizable 
operators. 

In the kincmatical region we are interested in, a more 
appropriate treatment of these contributions makes use 
of an expansion in powers of n ■ qA/q 2 ~ 0.37. This 
is similar to the approach adopted in Ref. |30l | for weak 
annihilation contributions to B — ► ir£ + £~. The SCETn 
operators obtained in this way are similar to those in 
Eq. (|3"0"|) . except that the soft operator is local. Keeping 
terms to second order in n-qh/q 2 , the spectator operator 
in this approach reads 

J s P = ^ E e q {^{qr la 4P L b) (31) 
+ ( V ~ (^«H DMPLb)} 

In this paper we will adopt the latter approach to the 
treatment of the spectator amplitude, working at leading 
order in n ■ qA/q 2 . We quote our results in terms of the 
first approach, which has an additional convolution over 
k-. However, it is straightforward to translate between 
the two approaches, simply by replacing 

J sp (fc_)~^ (32) 

<r 

below (in, e.g., Eqs. I|53|l and l|69p). The advantage of 
this approach is that the predictions are independent on 
the details of the matrix elements of the nonlocal soft 
operator (the B meson light-cone wave functions), but 
can be computed exactly in the soft pion limit to second 
order in the n ■ qA/q 2 expansion. 

In order to have a clean power counting of the transi- 
tion amplitudes, we divide the phase space of the B — > 
Ktt£ + £~ decay into several regions, shown in Fig. 

I) the region with one soft pion and one energetic kaon 
B -> K n Tr s , ~ A,Ek ~ Q and M 2 KlT - AQ. This 
region will be the main interest of our paper. 

II) the region B — > (K n ir n )K» describing decays into an 
energetic Ktt pair with a small invariant mass Mktt ~ A. 
This is dominated by one-body decays into a collinear 
meson B — > K*, followed by if* — > K n ir n . This region 
will be treated essentially the same way as a one-body 
decay. 

III) the region with a soft kaon and an energetic pion 
Ek ~ A, En ~ Q. The decay amplitude in this region 
is suppressed by A/Q relative to that in the other two 
regions I, II, and will be neglected in the rest of the paper. 

We will use the SCET formalism described above to de- 
rive a factorization relation in the region (I). The matrix 
elements of the nonfactorizable operators are parameter- 
ized in terms of soft nonperturbative matrix elements, in 
analogy with the B — > M n transition. We define them 
as complex functions of the momenta of the final state 
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FIG. 2: The phase space of the decay B -> Kni+r at <j 2 = 4 
GeV 2 , in variables (Mk-w, E^). The 3 regions shown corre- 
spond to: (I) soft pion ~ A; the shaded region < 0.5 
GeV shows the region of applicability of chiral perturbation 
theory; (II) collinear pion and kaon E^ ~ Q, Er- > 1 GeV; 
(III) soft kaon E K < 1 GeV. 



hadrons, with mass dimension zero 

{M n M' s \0^e%\B) = tl MU '(E M ,p M ,) (33) 
(M n M' s \0' ni \B) = (* MM '(E M ,p M ,) 

The matrix elements of the factorizable and spectator- 
type operators are given again by convolutions as in 
Eq. Ij25(l . with a different soft matrix element 



(M n M' s \Ji\B(v)) ~ / dxdzdk + b (l \z)J j (x,z,k + ) 

x(M' s \q k+ T s b v \B{v))(M n \q n!Ul Tcqn,uM (34) 
dxb sp (x)(p M (x) / dk-J sp (k-)(0\qk_Tsb v \B) . 



These factorization relations contain several new 
hadronic nonperturbative matrix elements, which we de- 
fine next. The new B — > n soft matrix element is defined 
in terms of the soft operator 



£ e -i^ +S (A|)F„(A,0)7^PR6.(0)(35) 



appearing in the term in Eq. I|29|l proportional to bm(z), 
and in the spectator operator. In the latter, one has to 
take into account that the light-like separation between 
the fields is along the direction n^, rather than as in 
the factorizable operators </f fact - 

The matrix element of the operator defines a soft 
function S as 

(n + {p«)\0,(k + )\B\v)) = -{g^ te^) P :S(k + ,t 2 ,p+) 

(36) 

with t = tubV — p-n- For simplicity of notation, we 
will drop the kincmatical arguments of the soft function 



S(k+, t 2 , Q whenever no risk of confusion is possible, and 
show explicitly only its dependence on the integration 
variable fc+. The matrix elements of the spectator op- 
erator in Eq. (|30|l are obtained from Eq. (|36[) . with the 
replacements n <-» n and ej^ v — * — £^„. 

The function S(k + ) is the B physics analog of a gen- 
eralized parton distribution function (GPD), commonly 
encountered in nucleoli physics |34j . The support of this 
function is the range —n-p^ < k + < oo, and its phys- 
ical interpretation is different for positive and negative 
values of k + . For k + > (the resonance region) the soft 
function gives the amplitude of finding a ud pair in the 
B° meson, while for k + < (the transition region), the 
soft function gives the amplitude for the b — ► u transi- 
tion of the B° meson into a 7r + meson. The soft function 
S(k+) is continuous at the transition point k + = |34j . 
which is important for ensuring the convergence of the 
k + convolutions in the factorization relation Eq. I|4U|I . 

We recall here the main properties of the soft func- 
tion S(k + ), which were discussed in Ref. |3^,|3^|. Time 
invariance of the strong interactions constrains it to be 
real. Its zeroth moment with respect to k + is given by 



X) 1 

dk + S{k + ,t\p+) = —n-pf T {t 2 ) 

pi 4 



(37) 



with frit 2 ) the B — > ir tensor form factor defined as 

(ir(p')\qi<T^b\B(p)) = f T (t 2 )( P ^ - p v p^) (38) 

Its iV-th moments with respect to fc + are related in a 
similar way to B — > ir form factors of dimension 3 + N 
heavy-light currents of the form q(n • iD) N b v [3fj| . 

In the soft pion region, chiral symmetry can be used to 
relate S(k + , t 2 , Q in the region k + > to one of the B me- 
son light-cone wave functions 4>?_(k+) defined in Eq. H28[l . 
according to 



S(fc+,t 2 ,c) = 



gfBrriB 



1 



4/tt v ■ p.* + A 



(39) 



Here g is the BB*ir coupling appearing in the lea ding 
order heavy hadron chiral effective Lagrangian |37l l38t 
I39L liiil ] . No such constraint is obtained using chiral sym- 
metry for S(k + , t 2 , C) in the transition region ( k + < 0). 

Collecting all the contributions, the amplitude for 
B — > M n M' s + leptons is given by a sum of factorizable 
and nonfactorizable terms, corresponding to the matrix 
elements of the SCETn operators in Eq. (|16fl . This leads 
to a factorization relation for such processes, which can 
be written schematically as 



A(B -» [M n M' s ] +leptons) = Cl (n ■ p M )( 



BMW 



(40) 



dz J dxdk+bi(z)Jj(x, z,k+)<f>M(x)S(k+\pM' 
dxb sp (x)(j)M(x) / dk + J sp (k + )S(k + ;p M <) 



This factorization relation has several important proper- 
ties |23l. First, the nonfactorizable contributions to the 
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decay amplitudes of semileptonic and radiative decays 
satisfy symmetry relations following from the universality 
of the soft matrix element ^ K7r . They contribute only 
to the decays B — > [M n M' s )h=-i£ + £~ into final hadronic 
states with total helicity — 1. Second, the amplitude for 
+1 helicity is factorizable, and given by a convolution as 
seen in the second term of Eq. (|40|l . Finally, the factoriz- 
able terms contain a new source of strong phase's, arising 
from the region k + < where the jet function develops a 
nonzero absorbtive part. This represents a new, factor- 
izable, mechanism for generating final state rescattering. 

Treating the spectator amplitudes in an expansion in 
powers of n ■ qA/q 2 according to Eq. (|3*T1) , the soft matrix 
elements are given by B — ► it form factors of dimension-3 
and 4 local operators. The leading order term contains 
the form factors of the vector current 

(n(p„)\u^b\B(p))=f + (t 2 )(p + p n )^+f_(t 2 )(p-p„)^ 

The f±(t 2 ) form factors appear in the matrix element of 
Eq. I|31[) in the combination /+ — /_. In the hard photon 
approach, the last term of Eq. I|40|) has the form A sp ~ 
frit 2 ) Jq dxb sp (x)(j)M(x), which follows from making the 
substitution Eq. ffify in this relation, and using Eq. IpTTjl . 

At subleading order in n ■ qA/q 2 , the form factors of 
dimension-4 currents u^aiD^b are also needed. They can 
be computed in the soft pion limit using chiral perturba- 
tion theory methods as discussed in Ref. |4lj . 

In the following section we derive the detailed form of 
these factorization relations for the B — > Kir£ + £~ decays. 



We start by recalling the results for the one-body de- 
cays B — > K* n l + £~ . The factorization relations for this 
case are well-known 0, E d, EH |H |27f and are given 
by (with i = V,A) 



H+\b — > K*) 
H^(B — > K*) 







(44) 
(45) 



3/ dzbfl(z)Cff (z) 
>o 



Note that the leading order spectator operator does not 
contribute to the decay with a transverse vector meson 
in the final state. 

The function Cjf'( z ) 

appearing in the factorizable 
term is defined as a convolution of the jet function with 
the light-cone wave functions of the K* and B mesons 



m B 



dxdk + J^(x, z, k + )(j)] 3 {k + )(j)j ( , {x\AQ) 



Using the result for the jet function Eq. 1)22(1 at leading 
order in a s (/j, c ), the integrals can be performed explicitly, 
and the function (j±(z) is given by 



/■BK" 



Nr. 



(47) 



with the first inverse moment of the B wave function 

0?(fc+) 



A" 1 
A B + 



dk. 



(48) 



III. FACTORIZATION RELATIONS FOR 

B -» Ktt£+£- 



The decay amplitudes B — ► K n it£ + £ into an energetic 
kaon and one soft pion can be parameterized in terms of 6 
independent helicity amplitudes H^' A \b — > K n ir) with 
A = ±1,0. They are defined as the matrix elements of 
the two hadronic currents in Eq. (jBJ 



H^ A \B K n n) = a»*(KM A \B(v)) 



(42) 



Working at leading order in 1/mf,, the helicity ampli- 
tudes can be written as a sum of nonfactorizablc and fac- 
torizable terms, arising from the corresponding SCETn 
operators in Eq. (|16|l 



H 



(B^K n n)= J2 H^ A) ' l (B^K 



(43) 



i— nf ,f ,sp 



The three contributions to each helicity amplitudes are 
computed as described in Sec. II. The nonfactorizablc 
terms are given in terms of the soft functions Cf KlT de- 
fined in Eq. (|33|l . and the factorizable and spectator con- 
tributions are given by factorization relations of the form 
shown in Eq. H4U|) . In this section we present explicit re- 
sults for the transverse helicity amplitudes. 



The corresponding amplitudes for the charge conjugate 
mode B — > K*£ + £~ are obtained from this by exchanging 
H+ <-*■ H-. The vanishing of the right-handed helicity 
amplitude at leading order in A/mb is a general result for 
the soft (nonfactorizable) component of the form factors 
in B — > M n , combined with the absence of the factor- 
izable contribution for this particular transition. This 
result is usually expressed as two exact symmetry rela- 
tions among the tensor and vector B — > V form factors 
at large recoil 0, . 

We proceed next to discuss the multibody decays B — * 
Kir£ + £~, in the kinematical region with q 2 ~ 4 GcV 2 . 
According to the discussion of Sec. II. A, the form of the 
factorization relation is different in the three regions of 
the Dalitz plot shown in Fig. [2J Our main interest is in 
the region I, with one energetic kaon, and a soft pion. In 
this paper we prove a new factorization relation for the 
transverse helicity amplitudes in this region. 

We consider for definitencss the mode B° — > 
K~ir + £ + £~ . Collecting the partial results in Sec. II. A, 
we find the following results for the transverse helicity 
amplitudes in this mode (with i = V,A), valid in the 
region I 



H^(B° - K-n+i+r) = H$ + &i,v\Hg (49) 



H®(B° K-TT+£+t 



(<) 



(50) 
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TABLE I: The transverse helicity amplitudes H± with i = 
V, A, for the different charge states in B — > K n ir£ + £~ decays, 



at leading order in A/Q. The building blocks H^ f ' , H\ ■ . 
are given in Eqs. 15H - 153j . 



H, 









B° -> ^-7r + ^+£- 
K s n°£ + r 




1 rrW 


B~ -> A'- 7 r°£ + ^- 




J_rr(i) 
J-ffW 



where the three terms correspond to the nonfactorizable, 
spectator and factorizablc terms in Eq. (|43|) , respectively. 
They are given by 



H, 



f ~ 2 
1 



Jk(£*+ - Pit) 



(51) 

(52) 



x / dzdxbiftiz) I dk+J\\{x, z, k + )S(k + )<f>K(x) , 
Jo '-pi 

H$ = ^f K (n-p K )(e*_-p«) (53) 
x / dxb$(x)<f>K{x) / dfc_J sp (fc_)5(fc_). 



The nonfactorizable operators contribute only to the 
left-handed helicity amplitudes, and are given by the soft 
functions . They are the same for both i = V,A 

amplitudes. Furthermore, the same soft functions would 
appear also in factorization relations for semileptonic de- 
cays into multibody states, such as B — > ir n ir£D. This 
universality is the analog of the form factor relations for 
the nonfactorizable amplitudes 0, @, well-known from 
one-body decays, to the multibody case. 

The factorizable operators give nonvanishing contribu- 
tions Hf to the right-handed helicity amplitudes. The 
appearance of these contributions is a new effect, specific 
to the multibody decays [23L l27j . On the other hand, the 
spectator operator contributes only to the left-handed 
helicity amplitudes. 

The helicity amplitudes for all other B — > K n ir£ + £~ 
decays can be obtained in a similar way. The results are 
tabulated in Table [I] 

The structure of these results displays universality of 
hard-collinear effects. This is manifested as the fact that 
all factorizable terms depend on the same x, k + con- 
volution 



dx 



dk + J\\ (x, z, k + )S(k + )4> K (x)(54) 



and all spectator contributions depend on the same fc_ 
integral 



I sp (n •<?)=/ 
J — i 



dfc_J sp (fc_)S*(fc_) 



(55) 



This universality is similar to that appearing in other fac- 
torization relations in exclusive decays. Examples are the 
relation between the rare leptonic decays B R — > £ + £~"f 
and the radiative leptonic decay B — > [2!j, and the 
relation among factorizable contributions in heavy-light 
form factors at large recoil, and the nonlcptonic B decays 
into two light mesons |3lj| . The integral I sp (n ■ q) appears 
also in the leading order factorization relation for exclu- 
sive semileptonic radiative decay B — > ir^lv |36l |. and 
could be determined from measurements of this decay. 
Treating the spectator amplitude using the n ■ qA/q 2 

expansion, the amplitude i?sp can be obtained as ex- 
plained from Eq. I|53|) , at leading order in this expansion, 
by the substitution Eq. 



4tt 2 



q^n ■ q 



-fKmBn-p K {£*--Pn)fT(t 2 ) (56) 



x / dzbil\z)^ K (z) 
Jo 

In the numerical estimates of this paper, we will use the 
leading order chiral perturbation theory result Eq. 1391) 
for the soft function S(k + ) in the resonance region k + > 



0. The contribution to the k + convolutions in H { 



(0 



and 



.Hsp from the transition region — n ■ p^ < k + < will be 
neglected, which can be expected to be a good approx- 
imation for very soft pions (the region < 500 McV, 
corresponding to the lower shaded region in Fig. |2J • We 
emphasize that, although the use of the chiral perturba- 
tion theory result for S(k+) is restricted to part of the 
region (I), the factorization relations proved in this pa- 
per are valid over the entire region (I) . Their unrestricted 
application requires a model for the soft function S(k+) 
whose validity extends beyond the limitations of chiral 
perturbation theory. 

With these approximations, additional universality 
emerges, connecting the amplitudes in this problem to 
other B decays, to all orders in the perturbative expan- 
sion at the hard-collincar scale. The factorizable helicity 

(i) 

amplitudes iJ f take a simpler form, and can be written 
as 



H 



1 



S R (p n ) I dzb<&{z)g K (z) (57) 



where the nonperturbative dynamics is contained in the 
factorizablc convolution defined as 



/-BK 



JbIk 
m B 



1 POO 

dx / dk+J\\ (x, z, k + )(f)^(k + )(j) K (xX58) 
o Jo 



The same function appears also in the factorizable con- 
tribution to the B — > K form factors at large recoil [Tfij . 
and in the factorization relation for nonlcptonic decays 
B -> KM with M = ir,K, ■ ■ ■ a light meson [3l|- The 
pion momentum dependence in Eq. 1)57(1 is contained in 
the function SniPn) given by [23| 



Sr{Pk) 



9 



Pr, 



Uvp w + A-iT B ./2 



(59) 
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with A = m B > - m B ~ 50 MeV. 

A similar result is obtained for the spectator amplitude 
at leading order in chiral perturbation theory, for which 
we find (treating the q 2 as a hard-collinear scale) 



4tt 



H $ = -^rfBfKm B (n-pK) 



g £--Pn 



(60) 



<t ' ' V U V ■ Pit + A 

pi poo 

x / dxb^{x)4> K {x) / dfc_J. p (ft_)0f 
Jo JO 

The k- convolution in this relation is identical to that 
appearing in the leading order factorization relation for 
radiative scmilcptonic decays B — > ^tv. 

Adopting the approach of expanding in n ■ qA/q 2 , the 
result for H^p requires the B — > ir form factors, for which 
one finds at leading order in HHChPT (H |H EJ 

fsruB 1 



m B Mt t ) = f+(t i )-f-(t t ) = -gi 



fir En 



(61) 



We will use these expressions together with Eq. 15611 in 
the numerical evaluations of Sec. V. 



IV. DECAY RATES AND THE FB 
ASYMMETRY 



The differential decay rate for B 
by (see, e.g. H2) 



Kirt+t 



1 



d 2 T 



r dq 2 dcos9 + dM 2 (yT dE 7r 2(4ir) 3 m 2 B m 5 b 



is given 



(62) 



x<^ 2 sin 2 + (I -Hi 



V\2 



\H Q A \ 2 ) 



A\2 



{ 

{l + cos 2 6 + )(\Hl\ 2 
+ 4cos6 + Rc(HYH A * - hXH^*)} 



H 



V\2 



\H 



A\2\ 



with T = G|a 2 /(327T 4 )|A t (s) | 2 mf. We denoted 9+ the 
angle between the direction of the positron momentum 
and the decay axis in the rest frame of the lepton pair, 
for a fixed configuration of the hadronic state Kit defined 
by (M Kn ,E„). 

Integrating over cos# + one finds for the forward- 
backward asymmetry (FBA) defined as in Eq. (JTJ 



Apb oc Re (HYH a * - HYHi 



(63) 



This defines a triply differential asymmetry depending 
on (q 2 , Mktt, E n ). Integrating also over E„ gives a dou- 
bly differential Afb depending only on (q 2 , Mk-k)- We 
denote them with the same symbol, and distinguish be- 
tween them by their arguments. 

The condition for a zero of the FBA can be writ- 
ten down straightforwardly using the expressions for the 
hclicity amplitudes in factorization given previously in 
Sec. III. The equation for the zero is different in the two 
regions (I) and (II) , according to the different form of the 



factorization relations in each of them. It is convenient 
to write this equation in a common form in both regions, 
as 



Re 



(c { p{M K „,q 2 )-a) 



0. 



(64) 



The quantity a stands for the contribution of the factor- 
izable and spectator type amplitudes, and in general is a 
function of all kincmatical variables (Mk„, q 2 , E n ). 

In the region (II) with a collinear kaon and pion, this 
correction is given by the (complex) quantity 



E K * n ■ q 



BK*eS 
J± 



/-BK* 



(65) 



where the factorizablc coefficient, (ff cff , which has im- 
plicit dependence on (MkttiQ )■> is defined by 



dzb&wff » = cfc?r* (66) 



The zero of the FBA in region (II) was considered previ- 
ously in Refs. 0, SIS El El > treating the problem as a 
one-body decay B — > K*£ + £~. 

In the region (I) with a soft pion and a collinear kaon, 
this correction contains two terms, arising from the spec- 
tator and the factorizable contributions, respectively We 
adopt everywhere in the following the leading order chi- 
ral perturbation theory results for the amplitudes given 
in Eqs. i|57|l and i|60|) . Working at tree level in matching 
at the scale fi = Q, but to all orders in the hard-collinear 
scale, one finds 



ai 



BKtt 



(67) 



(^Y\Sh(p,)\ 2 C, 



fBKsBKeff 
effSj 



\AE K J 



ICf 



BK-n\2 



The first term contains the contribution of the spectator 
amplitude, and depends on the charge states of the final 
and initial state through the superscript q = u,d, denot- 
ing the flavor of the quark attaching to the photon. The 
quantity /lip , given by 

a 2 r 1 

h$ = -^rfBfKm B (n-p K ) / dxbifix^Kix) (68) 
Q Jo 



X / dk-J sp (k-)(j>Z(k-) 

Jo 

depends on kinematic variables only through the explicit 
factor oin-px/q 2 , to the order we are working. For com- 
pleteness, we quote also the expression for this amplitude 
at leading order in n ■ qA/q 2 , which is used in the actual 
numerical computation of Sec. V 



4m 3 



q 2 n ■ q 



fBfKm B (n-pK) 



(69) 
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The first term in Eq. (|67|l contributes to Rc(ai) with an 
undetermined sign, depending on the unknown £j 3A_7r . 

The second term in Eq. H67(l is due to the right-handed 
hclicity amplitudes H + ■ Its dependence on is ex- 
plicit in SrIptt), and the remaining factors depend only 
on (Mkk, q 2 )- The factorizable coefficients in the numer- 
ator are defined as 



dzb<${z)Cf K (z)- 
dzb^(zKf K (z) 



1 



n -pk 
1 



CfCf Kofi (70) 



"CioCj 



BK 



(71) 



n -p K 

where we kept only the tree level matching result for 
^ir( z )- Numerical evaluation of these factorizable am- 
plitudes in the next section shows that the contribution 
of this term to Re(ai) is positive. 

To explore the implications of these results, let us as- 
sume as a starting point that the helicity amplitudes 
H± in both regions are dominated by the nonfactoriz- 
able contributions, proportional to ^ K7r (hi region (I)), 
and £± K (in region (II)). This corresponds to taking 
a = in both regions. In this approximation, the condi- 
tion for the zero of the FBA reads simply Re(ci ) = 0, 
which can be solved exactly. For the B — > M n transition, 
this condition reproduces the well-known result follo wing 
from the large energy form factor relations 

Since the zero of the FBA is related to the vanish- 
ing of the Wilson coefficient Re (c^ ), such a zero must 
be present also for B decays into multibody states con- 
taining one energetic kaon p3(|. In particular, the FBA 
in B — ► K n ir£ + £~ must have a zero at a certain point 
9o = Qo(Mkw) which depends only on the invariant mass 
of the hadronic system. Adding the second term in 
Eq. I|t)4|) shifts the position of this zero, and introduces a 
dependence on the pion energy, q^ = q^Mx-i^E^). This 
extends the well-known result for the zero of the FBA in 
B — > K*£ + £~ to multibody hadronic states. 

It is interesting to comment on existing computa- 
tions of the decay am plit udes and FB asymmetry in 
B -> Kn£+e- HlaElllil, which keep only the K* res- 
onant amplitude. Of course, this is justified in the region 
(II), where the pion and the kaon are collinear. However, 
in the region (I) this contribution is in fact paramctrically 
suppressed, since by leading order soft-collinear factoriza- 
tion, the K*K n TTs vertex does not exist at leading order 
in 0(A/m b ). 

Computing the factorizable corrections in the region 
(I) parameterized by the term a\ requires that we know 
the nonfactorizable soft function C^ Kv . There are sev- 
eral possible ways of determining ^ Kt! from data. For 
example, according to Eq. the helicity amplitude 

receives no factorizable or spectator contributions. 
Assuming that it can be isolated, its measurement would 
give a clean determination of ICf^l- Another method 
involves measuring Hzf^ in decays with a neutral kaon 
in the final state, for which the spectator contribution 
is small (see Table I). We assume that the factorizable 
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FIG. 3: Plot of the zero in the Forward-Backward Asymme- 
try as a function of the parameter F of Eq. 17311 . for several 
values of M Kl . The top (green) curve has lowest M Kl and 
the bottom (blue) one has highest Mktv- Fq corresponds to 
the value of F in the absence of factorizable and spectator 
corrections. The maxima of the curves, labelled by F„, lie to 
the left of Fq for all relevant values of Mjf,. 



coefficients can be computed in perturbation theory. We 
postpone a detailed numerical analysis for Sec. V, and 
discuss in the following general properties of the zero of 
the FBA which are independent on the details of the 
hadronic parameters. 



A. Qualitative discussion of the zero of the FBA 

Before proceeding with the details of the numerical 
study, we would like to discuss some of the qualitative 
properties of the zero of the FBA in the multibody decay 
B — > Kit£ + £~. The general behaviour of the solution 
can be seen by studying the solutions of the simplified 
equation 



C 9 + 2m b ^-C 7 - a(Mj(ir) = 



(72) 



where a{Mx-n) denotes the factorizable correction. In 
this simplified version of Eq. H(j4|) we have neglected ad- 
ditional dependence of a on q 2 , which is adequate if one 
is interested in the qualitative change in the zero in the 
FBA at fixed q 2 . We also have neglected here the radia- 
tive corrections. They introduce small logarithmic de- 
pendence on Mk-ki and do not change the qualitative 
features of the solution. 

It is convenient to write the equation Eq. I|72l) in an 
equivalent form 



n • q 

2 

q- 



2m b G 



(C 9 -a(M & ))E*^ (73) 



m B 



The solution to this equation gives the location, q 2 ,, of 
the zero of the FBA: 



M 2 
F-l 



(74) 
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The condition that the FBA-zero lies in the physical re- 
gion, < <7q < {m B — Mk-r) 1 , imposes constraints on 
F = F(M Kn ). The upper bound g§ < (m B - M K7T ) 2 
implies that F(Mk-r) > 1, while from < g§ we learn 
that 



F(M K „) > 1 



> 1 



(75) 



In terms of the correction a the condition for the exis- 
tence of a FBA-zero is therefore 



a{M K „) < C 9 



2mbmB 

? 2 _ ]\J2 



-c 7 . 



(76) 



If the function F = F(Mk-r) is roughly constant and 
satisfies the condition in 1|75[) . then Eq. I|74|) gives that the 
zero of the FBA decreases with increasing Mkr- These 
conditions hold if cl{Mk-k) "C Cg. Since we expect the 
correction term to be small we also expect the FBA-zero 
to decrease as Mk-r increases. 

We can gain further insights into the solution by con- 
sidering (/q as a function of F for fixed Mk-r in Eq. (|74|l . 
Fig. shows plots of q 2 vs. F for several fixed values 
of Mk-r- The sequence of curves moves down with in- 
creasing Mktt, which just restates the observation that 
at fixed F the zero decreases with increasing Mk-k- The 
point F ~ Fq corresponds to a = 0. The maxima of 
the curves are at F^(Mk-r) = (1 — M k„ / 'raj) -1 , and, for 
physical values they lie to the left of i*o, that is, F* < Fq. 
To see how q 2 depends on Mk-r when F = F(Mkr) is 
not constant, consider as starting value a point on the top 
curve. An increase in Mk-r first moves the point down to 
a lower curve (as if F were constant), and then also along 
the lower curve to a different value of F. The region to 
the right of F* is most interesting since we expect the 
physical function to lie in a region of F close to Fq. In 
this region, if F increases with Mk-r then q 2 , decreases 
with (increasing) Mk-r- The opposite is not necesarily 
true: whether q^ increases with Mk-r or not depends on 
how steeply F decreases with Mk-r- 

Next, we would like to understand what the effect of 
changing a — > a + Sa is, corresponding to adding correc- 
tion terms sequentially. A shift SF > for fixed Mk-r 
corresponds to moving to the right along a fixed curve. 
In the region to the right of F* this decreases q 2 . Note 
that 6a = (2m b C 7 /m B C 9 )SF and C 7 /C 9 < 0. Therefore, 
an increase in a gives an increase in q 2 ,. 



V. NUMERICAL STUDY 

We investigate in this section the numerical effects 
of the new right-handed amplitude on the position of 
the zero of the FB asymmetry. As explained, we treat 
separately the decay amplitudes in the two regions (I) 
and (II), and ignore the contribution from region (III). 
For definitcness, we consider here the mode B° — > 
K~ir + £ + £~ for which the soft pion detection efficiency 
is better than for the neutral pion modes. 
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FIG. 4: Plot of the position of the zero of the forward- 
backward asymmetry qo = q$(Mx ) as a function of the in- 
variant mass of the Ktt system, obtained by neglecting the 
factorizable contributions to the helicity amplitudes, for dif- 
ferent values of the renomalization point fi. 



The decay amplitudes in the collinear region (II) will 
be represented by a Brcit-Wigncr model, as 



hY{b° -> K-ir+i+r) = h v _{B -> K*) 

X9k*kA £ - ■ Pr)BWk'(Mk^) 
H^(B° -> K-ir+e+r) = h A {B -> K*) 

XgK*KTr(e~ ■ Pr)BWk>(M K k) 



(77) 
(78) 



and H+' A = 0, with li^ A the one-body helicity ampli- 
tudes for B — > K*£ + £~ given above in Eq. The 
Breit-Wigner function corresponding to a K* resonance 
is defined as 



BWk-(M) 



1 



M 2 



M 2 



iMk*Tk* 



(79) 



Finally, the K*°K + ir coupling with a charged pion can 
be determined from the total K* — ► Ktt width, T = 
g K » Kn p%/(l&'Km K ,), with the result g K -K-R = 9.1. 

The factorization relations in region (I) require the 
nonfactorizablc amplitude . In the absence of ex- 

perimental information about this quantity, we adopt a 
K* resonance model for it, defined as 



£ K *{M Kn ,E„) 



(80) 



*9k*kA £ - ■Pr)BWk'(Mkr) 



where the kinematical factor n-pK* = 2Ek* = ^"( m s — 
q 2 ) can be chosen corresponding to an on-shcll K* me- 
son. For the kinematical dependence of the soft function 
f (q 2 ) we adopt a modified pole shape [THL lisj 



cr (o) 



l-1.55^ + 0.575(^-) 2 



(81) 



and quote results corresponding to the two values 
Cf K *(0) = 0.3 and 0.1. These two choices should cover 
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both cases of soft-dominated, and hard-dominated tensor 
form factor. 

We will use this model to define a FBA differential in 
(g 2 , Mk-k), integrated over the pion energy E„. Separat- 
ing the contributions from the regions (I) and (II), this 
is given by 



can be expressed in terms of three phase space integrals 
-^0 1,2) arising from region (I), and another integral I in 
region (II), defined as 



A FB (q 2 ,M Kv ) = / dE„Re [HYH** - H^Hf] 
J(i) 

+ [ dE^Re [H^H M ] (82) 
J(ii) 

The integration over E w can be simplified by approxi- 
mating n ■ px f» niB — n ■ q in region I. Then the result 

I 

The zero of the FBA is given by the solution of the equation 

Re [c^ y) - a sp - Of] 



i 



dE* !l + ' P l l l ■ 3 = 0,1, 2 (83) 



dE 7r \e+ -Pir\ 



(84) 



E out 



These integrals depend implicitly on (Mk„, q 2 ). Numer- 
ically, for E^ ut = 500 MeV we find I = 0.019 GeV 3 , 
h = 0.047 GeV 2 , h = 0.12 GeV and I = 0.12 GeV 3 , at 
M K7r = 1 GeV and q 2 = 4 GeV 2 . 



0. 



(85) 



where a sp and af denote the factorizable contributions, arising from the spectator amplitude, and from the factorizable 
(both in one-body and in the two-body amplitudes). The a sp coefficient is given by 



9 



3 U n ■ Vk*9k* Cl h + 



-[M Kw -M K *+iM K *T K .] 



(86) 



where hip is defined in Eq. I|69[) . 

The factorizable term contains contributions from the one-body decay amplitude, and from the new right-handed 
amplitude appearing in the two-body mode 



Of 



2m% 



/■BK* cS 



2 4 
9 ™>B 



n-p K ,n-q Cf K ' I + I Af 2 g 2 K , K Jn ■ p K ) 2 (n ■ p K 



C 



/■BK/-BKc« 



(M K „-M K ,f 



h 



(87) 



We compute the factorizable matrix elements using the 
leading order jet functions from Eq. 1)22(1 . This gives for 
the integrals of the factorizable functions (j(z) 



t BK 
/-BK* 



and for the effective ones 



■na s C F Sb!k 1 f 1 J z 4n 
N a rns A B + Jo 1 

^ i tidz^ 



4>k{z) 



TTO. s Cf IbI 



/-BKcft _ /-BK 



cir 



/■BK'cti 



/■BK* 



2m b ((i)n ■ q Cf 

q 2 Cf 
e-uE K C 2 j BJf 
?7i fc n • g C§ ff J 
2e u E K , C 2 BK * 
1 + ^n~Tcf Ij± 



where 



/f * = /o 1 dzt x (2, m c )0^) / Jo , 
fZdzt ± (z t m c )<l%.(z)/fZdz&¥ 



jBK 
'J_L 



(88) 
(89) 



(90) 
(91) 

(92) 
(93) 



For the computation of the integrals we use the 
light-cone wave functions 



(j) K {x) = 6xxfl + 3a 1K (2x - 1) + ^a 2K [5(2x - l) 2 - 1]) 



(94) 



and analogous for <fi K ,(x), with coefficients aj~ K ». The 
values of the first two Gcgcnbauer moments are given in 
Table [H] Also listed there are the remaining hadronic 
parameters used in the computation. 

The resulting values of the factorizable matrix ele- 
ments are tabulated in Tablc lTTll which also lists the effec- 
tive Wilson coefficients Cfl (computed at \i = 4.8 GeV, 
g 2 = 4 GeV 2 ). The effective matrix elements (f KeS 
and (ff cff depend on (Mx-mq 2 ), partly through im- 
plicit dependence of the integrals lf K and iff . To 
gain some understanding of the relative importance of the 
terms that contribute to the effective matrix elements, 
we quote these integrals at Mk-k = mx* , q 2 = 4.0 GeV 2 : 
lf K = -0.704 - 2.564i and Iff* = -0.566 - 2.67i. 

It is straightforward to estimate the numerical impor- 
tance of each term in the correction a appearing in the 
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Cf *(0) . 
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FIG. 5: Plot of the position of the zero of the forward-backward asymmetry go — Qoi^Kn) as a function of the invariant mass 
of the Kn system. The plots show the change in the position of the zero due to the spectator and factorizable amplitudes, 
Eq. ED, f°r three values of the pion energy cut-off E^ ut = 300, 500 MeV and E^ ut = 700 MeV separating regions I and II. The 
dotted (blue) line denotes the position of the zero in the absence of the factorizable and spectator contributions (for /i = 4.8 
GeV). The non-factorizable matrix element is taken to be £f K " (0) = 0.3 (left) and C± K * (0) = 0.1 (right). 



TABLE II: Input parameters used in the numerical compu- 
tation. 



mf 


4.68 ± 0.03 GeV [43] 


Ik 


160 MeV 


fh c (fh c ) 


1224 ± 57 MeV [44] 


fl* 


175 MeV 


a s (Mz) 


0.119 


aiK 


0.3 


As+ 


350 MeV 


aiK 


0.1 


9 


0.5 [45] 


a{ K , 


0.2 


Sb 


200 MeV 


a 1K* 


0.1 




-0.0106 + 0.0174i [4£1 




9.1 



equation Eq. (|64|) for the zero of the FBA. To this end 
we evaluate at Mj(w = m K',(l 2 = 4.0 GeV 2 , and use 
££ ut = 500 MeV and K ' (0) = 0.3 in the model of 
Eq. JHU- We obtain (all in GeV units): 

a sp (0.016 + 0.0Ui)(M Kn - M K , + iM K *T K -) , 

Of ~ (0.533 + 0.32H) (95) 
+ (0.004 - 0.003i)(M Kn - m|.) 2 . 

We have kept explicit the rapidly varying dependence on 
the inverse Breit-Wigner function, so we may get some 
idea of the relative size of the coefficients. The first term 
in the factorizable amplitude can be regarded as a nega- 
tive correction to Re(Cg ff ) in of about ~ 10%. Thus 
it effectively shifts the zero of the FBA upwards by the 
same amount. The second term in at and a sp are negli- 
gible on resonance but may be important at large Mk-k- 
Using this model, we solve for the zeros of the FBA in 
B — > Kn£ + £~ decays, finding q$ for given [Mx-n, E^ nt ). 
The results are shown in Figs.^HSJ In Fig. 0] we plot the 
result for q^Mx*) obtained by neglecting the factoriz- 
able and spectator terms (the solution to Re (c^ ) = 



0) for three values of the renormalization scale, fi — 
2.4, 4.8, 9.6 GeV. We used here NNLL results for the Wil- 
son coefficients and the 2-loop matrix elements of the 
operators obtained in Ref. ^tJ- The position of the 
zero at threshold is 

($\m k «=m k+ m„ = 3.75^-i 2 GeV 2 (96) 

where the uncertainty includes only the scale depen- 
dence. This result depends only mildly on Mk-k, as seen 
from Fig. 0] 

In Fig.[3]we show also the effect of including the factor- 
izable and spectator terms in Eq. I|85|l , for three values of 
the cut-off on the pion energy E^ ut = 300 MeV, 500 MeV 
and 700 MeV separating regions I and II. The parame- 
ters used in evaluating this plot are listed in Tables [Q] 
and II I II For comparison, we present in Fig. results 
with Cf K *(0) = 0.3 (left) and (f K ' (0) = 0.1 (right). 
The latter choice effectively amplifies the factorizable and 
spectator corrections, and can be taken as a conservative 
upper bound of these effects. The dependence of the re- 
sults on E^ ut is an artifact of the separation of regions 
discussed above and indicates the uncerainty in this pro- 
cedure. For this reason we have taken rather extreme 
values of E% ut . The overall trend of g 2 (A/x7r) decreasing 
towards the right of the plots is readily understood from 
the qualitative discussion in Sec. IIV Al it follows from 
the correction term a being small. Similarly, that the 
inclusion of spectator and factorizable corrections tends 
to increase the value of q$ for fixed Mk-k follows from 
positivity of a. 

The results show a marked dependence (especially for 
small Mktt) of the zero position on the pion energy cut- 
off E^ ut which separates the regions (I) and (II). This 
is essentially due to the dominance of the factorizable 
contribution in region (II). A conservative way to use 
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our results is to take the smaller value of E^ ut = 300 
MeV, for which the chiral perturbation theory result can 
be expected to be the most precise. 



TABLE III: Results for the effective Wilson coefficients and 
factorizable and spectator matrix elements. The values of the 
effective Wilson coefficients are at the scale fi = 4.8 GeV and 
q 2 = 4 GeV 2 . The factorizable matrix elements are computed 



at the scale fi c = 1.5 GeV. 



[Cg ]nnll 4.579 + 0.082i 
[C? ff ] N NLL -0.388 - 0.020i 



0.036 
0.035 

(0.135 + 0.124i)GeV 3 



So far our considerations were restricted to the case of 
B° — ► K~ir + £ + £~ decays. Going over to the CP con- 
jugate mode B — ► K + ir~£ + £~ , the position of the zero 
could change because of direct CP violation present in the 
spectator amplitude , which is furthermore enhanced 
by a 47T 2 factor (see Eq. JHsUO- Denoting the factorizable 
corrections analogous to a sp and df for the CP conjugate 
mode with a,, we find 



■ -(0.016 - 0.014i)(MjL.-M£. +iM K .T K ») , (97) 



such that the CP asymmetry in the position of the zero 
is induced through the finite K* width, and is small. 
Beyond tree level, such an effect will be introduced at 
order a s (Q) through matching corrections to b sp . 

We consider next another observable, the slope of the 
curve for the zero of the FBA in Fig. [5] From Eq. (|74|l . 
this is given by 



ML 

K 71 



F- 1 \(F- l) 2 



m%\ dF{M*Y 



VI. CONCLUSIONS 

We studied in this paper the helicity structure of the 
exclusive rare B — > Kir£ + £~ decays in the region of phase 
space with one energetic kaon and a soft pion. In this 
region the helicity amplitudes are given by new factor- 
ization relations, containing an universal soft matrix ele- 
ment, and a new nonperturbativc matrix clement for the 
B — > 7r transition analogous to the off-forward parton 
distribution functions. 

The most important difference with the B — > K*£ + £~ 
decays at large recoil is the appearance in the multibody 
case of a nonvanishing right-handed helicity amplitude 
B — > {Kn]i l=+ i£ + £~ at leading order in A/m&. This can 
be computed in factorization, in terms of the B — > it off- 
forward matrix element of a nonlocal hcavy-to-light op- 
erator. In the soft pion limit this nonperturbativc matrix 
element can be computed in chiral perturbation theory, 
and is related to the B meson light-cone wave function 

M 

We explored the implications of these results for the 
existence of a zero of the forward-backward asymme- 
try of the lepton momentum, pointing out two new re- 
sults. First, the FBA has a zero also for nonresonant 
B — > Kir£ + £~ decays, occuring at a determined value of 
the dilepton invariant mass q^Mic-n), depending on the 
hadronic invariant mass Mjf T . Second, there are calcu- 
lable corrections to the position of the zero, which can 
be computed in factorization. We use the factorization 
relations derived in this paper to compute these correc- 
tion terms. We present explicit numerical results work- 
ing at leading order in chiral perturbation theory |23j . 
and show that the results for the zero of the FBA in 
B — > [Kit]k*£ + £~ hold to a good precision also in the 
nonresonant region. 



F(M Kn ) is defined in Eq. l(73)) and depends on the Wilson 
coefficients C^g, also on the factorizable contributions 
a(MKn)- The last term contributes through the Mk-k 
dependence of a(MKn), and is 



m B 



2mijCj 



dM 2 Klr 

K 7T 



0.02, 



(99) 



where we used the result Eq. for a sp (Afx7r) and ne- 
glected the tiny contribution from a/. The contribution 
of this term to Eq. I|98(l is multiplied with a factor of or- 
der 1-2. Thus, even assigning this estimate a conservative 
error of ~ 200%, its contribution to the slope Eq. I|98(l 
for values Mk-k ~ 1 GeV, is negligible compared to the 
first term depending only on F (recall that F ~ 4) . This 
is also seen in the curves in Fig. \5\ whose slopes are es- 
sentially the same for all choices of the hadronic param- 
eters considered. This shows that a measurement of the 
slope of the zero could provide a useful source of infor- 
mation about the Wilson coefficients 6*7,9, but without 
the hadronic uncertainties associated with the absolute 
position of the zero. 



APPENDIX: EFFECTIVE WILSON 
COEFFICIENTS 



We collect here for convenience the expressions for the 
effective Wilson coefficients used in the numerical study 
of Section IV. Working to NNLL order, they are given by 



-C 9 - (Ci + y )(8G(m c ) + 1) 
-C 3 (8G(m c ) - 1g(0) - yG(m h ) + A) 
+G 4 (4G(0) - *-G{m c ) + yG(m h ) + 
-G 5 (8G(m c )-4G(m 6 )-^) 
-C G {\G{m c )-- i G{m b )+ 2 -) 



(A.l) 



7-,(9)/ 9\ 



•£[2Ci(JTV) ■ 
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and 



Cf = C 7 - ^0 3 - ^0 4 + ^C 5 + ^0 6 (A.2) 



-gPa^VH^'V)] 

They are expressed in terms of the modified Wilson coef- 
ficients Ci_6, which are defined by expressing the oper- 
ators 0i_6 of Ref. |24| in terms of the basis of 0] using 
4-dimensional Fierz identities. They are given by Q 





=\* 




02 — 02 


c 3 


= 03- 




+ 1605 " 


c 5 


= 03- 




+ 40 5 - 



(A.3) 



The functions F^fa 2 ), F^\q 2 ) appearing in the 2-loop 
matching conditions are listed in Eqs. (54)- (56) of the 

second reference in Ref. j^. The functions Fg 7 ' 9 \q 2 ) 
are given in Eqs. (82), (83) of Ref. Q. 

We list here the functions f v and r appearing in the 
expression of the Wilson coefficient Eq. © 



-2 log 2 



1 i 2 m b 

o lo § -r 

2 [i z 
2Li 2 (l 



5 ni 
2 /j,- 



2 log -i log- 

a; . 3oj — 2to(, 

— ) + log — 

nib n^b — w wi& 



— + 6. 
12 



a s 0i 



4tt 



a; a; 
log — 

nib — w ^6 



log- 



(A.5) 
(A.6) 



where C,. are the Wilson coefficients in the operator basis 
of Ref. j24|. The 2 ; coefficients coincide with the Wil- 
son coefficients in the basis of Ref. ^(| , but are different 
beyond leading log approximation. The relation between 
the two sets of coefficients can be found in Refs. 
The effective Wilson coefficient 0g ff is given by 



Ccff r ' i s~i ^ s~i 

8 — C '8 + - ■ 

The one-loop function G(m q ) is given by 



(A.4) 



f 1 / — q 2 xi\ — x) + m„ — ie 
G(m q ) = J dxx{l - x) log (-1-^ -L q - 



In the numerical evaluation of the Wilson coefficient 



pole 



we replace the MS mass mb(fi) with the pole mass 
using the one-loop result 



rrib(n) = ml° lc { 1 



^(-61ogi--4) 
47r mb 



and keeping only the term linear in a s . 

Finally, we give here the function t±(x,m c ) appearing 
in the Wilson coefficient of the sublcading O(A) SCETi 
operators. This is given in Eq. (27) of Ref. 0, which we 
reproduce here for completeness 



t±(x,m c ) = -^Ji(m c ) + -^-^[Bo{xm 2 B + xq 2 7 m c ) - B (q 2 ,m c )} (A. 

with 



/ 4m 2 1 
B (s, m c ) = -2\ - 1 arctan — ^=^= (A.9) 

2m 2 c 
x{m 2 B — q 2 ) 1 



h(m c ) = 1 + = , 2 c 2 J L x (x+) + - L l {y+) - L x {y_)} 



The function Li(x) and its arguments are defined as ACKNOWLEDGMENTS 
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The existence of a zero of the FBA in inclusive B — > 
X s £ + £~ decays was noticed and studied in Refs. fiolllll . 
In this paper we address only the semi-inclusive decays 
B — > K n n with low hadronic invariant mass, for which 
the inclusive methods of Refs. [Tol . llll are not applicable. 



